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Abstract 

We generalise a recent derivation of the relativistic expressions for momentum 
and kinetic energy from the one-dimensional to the three-dimensional case. 
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1 Introduction 

We have recently shown pQ how to construct, given a velocity composition law, the expres- 
sions of kinetic energy, momentum, the Lagrangian and the Hamiltonian for a particle, 
in a general mechanical theory satisfying the principle of relativity, and in which elastic 
collisions between asymptotically free particles exist. For reasons explained in that paper, 
the discussion in reference was restricted to the case of one spatial dimension only. In 
the present article we extend the treatment to three dimensions. 

Let us briefly review the key elements upon which the derivation in reference 
is based. The starting point is the definition of kinetic energy for a particle, as a scalar 
quantity whose change equals the work done on the particle. Mathematically, this amounts 
to requiring the validity of the fundamental relation 

dT(u) = u • dp{u) , (1.1) 

where T(u) and p(u) are, respectively, the kinetic energy and the momentum of a particle 
with velocity u. For a system of non-interacting particles, kinetic energy is then necessar- 
ily additive, since work is, so one can easily write down a formula that expresses energy 
conservation during an elastic collision. On requiring that this holds in an arbitrary iner- 
tial frame, it follows that another quantity is conserved, in addition to energy — a vector 
one, that we identify with momentum. As we shall see, this quantity is linked to kinetic 
energy through a simple equation containing a matrix of functions tpij (the indices % and 
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j run from 1 to 3), which is uniquely determined by the composition law for velocities. 
This, together with equation (jl.lj) above, allows one to find the explicit dependence of 
both energy and momentum on velocity. 

We begin by collecting, in the next section, some results concerning the velocity com- 
position law in three dimensions. Section |3] contains the general outline of the method, 
which is then applied in section 0] to the cases of Newtonian and Einstein's dynamics. A 
short summary of the results, together with an assessment about their suitability from a 
pedagogical point of view, is given in section El The Appendix contains two mathematical 
relations that are not used in the body of the paper, but that might nevertheless turn out 
to be useful in further possible developments. 



2 Kinematical preliminaries 

Suppose that a particle moves with velocity u with respect to a reference frame /C. If K, 
moves with velocity v with respect to another reference frame /C, the particle velocity u 
with respect to K, is given by some composition law 



u = <&(u, v) 



(2.1) 



(It is important to appreciate that this relation contains vectors belonging to two different 
spaces. Not only are the basis used to write the vectors v and u, and the vector u, 
different; they even span different rest spaces, namely those of /C and JC, respectively.) 
The relativity principle requires that ()2.1|) give the composition law of a group, i.e., that: 



$(u,0) = $(0,w) = u , ; 

Vti, 3u such that $>(u, u) = &(u', u) = u ; 
<&(<I>(ii, v), w) = &(u, w)) , Wu,v,w . 
Note that, although in the Galilean case, 

<fr(it, v) = U + V , 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



the composition law is commutative, this is not the case in general, unless velocities are 
collinear. For example, the relativistic law [2] 



<fr(ii, v) 

where, as usual, 



1 + u • v /c 2 



u 



j(v) 



+ v 



U ■ V 



7(v) 



+ 1 



y/l ~ V 2 /C 2 



(2.6) 



(2.7) 



denotes the Lorentz factor, is not commutative. 

From the composition law (J2.1|) one can define a matrix whose components are, in 
Cartesian bases, 



<PijW 



(2i 



v = 
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Equation (J2.3j) then imposes that <fij{0) = 6^. In fact, for the Galilean composition law 
(12. 51) . we have (pij(u) = 5^ for all u. On the other hand, in the relativistic case (|2.6j) . 

( p ij (v) = S ij -^. (2.9) 

3 General analysis 

Let T(u) be the kinetic energy of a particle with velocity u in an inertial frame /C. During 
an elastic collision between two particles, energy conservation requires that 

71(1*0 + T 2 (u 2 ) = 71(w'0 + T 2 {u' 2 ) . (3.1) 

(Of course, the kinetic energy will also depend on the particle mass; we keep track of this 
dependence with the indices 1 and 2 on T.) 

With respect to another inertial frame 1C, in which K, moves with velocity v, the 
particle velocities are Ui = $(ui,v), u 2 = &(u 2 ,v), u[ = <f>(u[,v), and u 2 = &(u' 2 ,v). 
Conservation of energy in K, then implies 

71(w0 + T 2 {u 2 ) = IKu'O + T 2 {u' 2 ) , (3.2) 

where we have used the same functions 71 and T 2 as in (j3.1|) . because of the relativity 
principle. 

The expansion of T(u) around v = is 1 

T(u) = T{u) + vi Vij {u) + 0{v 2 ) . (3.3) 

Doing this for each term in equation (|3.2|) . then using equation ()3.1)1 and considering the 
terms of first order in v, we find the following additional conservation law: 

+ ^(-0 = Vy(«i) + ^K) • (3-4) 

Thus, one arrives at the conclusion that the vector p(u), whose components are 

Pi{u) = (pij[u) — , (3.5) 

is conserved during a collision, in addition to energy. In the one-dimensional case, this 
quantity can be identified with linear momentum |]Q, and we suggest doing the same in 
three dimensions. 2 Note that, with this definition, linear momentum turns out to be a 
one-form rather than a vector, which is very satisfactory from a formal point of view. 

^^We adopt the convention of summing over repeated indices. 
2 Actually, the most general conserved quantity has the form 

. , dT(u) , . 

Xij ipjk(u) —q^- h m T(u) + Vi , 

where Ay, fit and Vi are a tensor and two vectors that do not depend on u. The simplest choice A^ = dij, 
[ii = vi = corresponds to the requirement that space be isotropic, and is the one which leads to the 
standard expressions for T(u) and p(u). 
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If we know the function T(u), we can find p. If we do not already know T(u), we can 
define it by requiring that it satisfies the fundamental relation (jl.ljl . which expresses the 
equality between the work done on the particle and the variation of its kinetic energy |T] . 
On rewriting the differentials dT and dp in (jl.lj) in terms of dui one gets 

^W =M .%W. (3.6) 
dui 3 dui 

Taken together, equations (|3.5|) and (|3.6|) allow one to determine T(u). Using again (|3.5|) . 
one can then find p(u). 

The free particle Lagrangian must satisfy the relation 

PiW = ■ (3-7) 

Using equation ()3.5j) . we obtain 

dT(u) 

dL{u) = Vij (u) d Ui . (3.8) 

Obviously, it is only for (p^ = 5ij (i.e., in Newtonian dynamics) that L = T + const — a 
feature already emphasised in reference pQ. 

Turning now to the Hamiltonian, we need only notice that ()3.6|) gives, basically, half of 
Hamilton's equations of motion for a system with Hamiltonian H(p) = T(u(p)) + const. 
Indeed, 

^ = d Uj (p)dT(u) = dH(p) 

dpi duj dpi 

or, symbolically, u = V p iL This allows us to identify H(p) with T(u(p)), up to an 
additive p- independent constant. 



4 Isotropy 

The previous discussion was general, in the sense that it was based only on the principle 
of relativity and on the hypothesis of space homogeneity (implicit in our use of inertial 
systems). With the further requirement that space be isotropic, one can restrict (pij(u) 
to having the functional form 

<fij(u) = 5ij + f(u) Ui uj , (4.1) 

where / is an arbitrary function of the magnitude u of u. This follows immediately by 
considering that no other vector except u can be used in writing ipij. In fact, even the 
class f!4.1|) is too wide, because relativity, homogeneity, and isotropy together, force / to 
be a constant. 3 Simple physical considerations j3] then require that such a constant be 
non-positive, so we shall write from now on 

(Pij{u) = 5 i:j -Kui Uj , (4.2) 

3 Since / depends only on the magnitude u of the velocity, this result can be established simply by 
comparison with the one-dimensional case ^ |3] . 
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where K > 0. The cases K — and = 1/c 2 correspond to the Galilei and Einstein 
composition law. 

From equation ()3.6|) we find 

dT d(pjUj) . . 

Ui— = Ui ^ 31 - PiUi . (4.3) 

Inserting ()3.5|) with the form (|4.2|1 for into ()4.3|) . we obtain 

2 Mi — = Ui — [uj K u Uj 7^— + K u Ui——. (4.4) 

OUi OUi \ OUj OUj J OUi 



Using the mathematical identity 



d „ 2 d „ d 



where £ := lnu 2 , equation (|4.4jl can be rewritten as 

(2 + ^)^ = 2(1-^0. (4.6) 

This is a simple differential equation for T, that can be solved by first finding dT/d£: 
dT . /I /■ 2 + Ke£\ Ae« 

where A is an arbitrary constant. By a further integration we find, for K ^ 



T{u) = A,K m + B ; (U 

'1 — K u 2 ^ 



while, for K = 0, 

T{u)=Au 2 + B. (4.9) 

In both cases, B is a further integration constant. On requiring that T(0) = 0, we can 
write .B in terms of A, which can then be expressed using the more familiar parameter m 
- the particle mass. We thus obtain, for K ^ 0, 



m/K 
'1 -Ku< 



T{u)= - 1 m/K; (4.10) 



and, for K = 0, 

T(n) = ^mu 2 . (4.11) 

Note that (|4.11j) coincides with the limit of (|4.10|) for K — > 0. For if = 1/c 2 one recovers 
the standard expression for kinetic energy in Einstein's dynamics, 
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It is then a straightforward exercise to find the usual expressions for momentum, the 
Lagrangian, and the Hamiltonian, using equations (|3.5j) . (|3.8j) . and the remarks at the 
end of section El The results are the same as in the one- dimensional case, so we do not 
list them here. 

In closing this section, we note that the calculations can be simplified by arguing that 
T[u) can only be a function of u, and that p(u) must be of the form p(u) = a(u)u, with 
a an unspecified function of u. However, we have preferred not to rely on these results, 
because they need not be derived separately, but are actually consequences of equations 
(P , lEOfll . and (P . 



5 Conclusions 

In this paper we have extended the derivation of relativistic energy and momentum given 
in pQ , from one to three dimensions. Although nothing changes in the results, it is obvious 
that the discussion is not as elementary as for the one- dimensional case, since it requires 
some familiarity with multivariable calculus. Hence, this material is not suitable for an 
introductory course, contrary to what happens for the treatment in reference pQ. It can, 
however, be presented in a standard course on analytical mechanics, at the level, e.g., of 
reference |3J. Indeed, we believe that students would benefit from being exposed to this 
approach, which relies on a very small number of physical hypotheses, and has therefore 
the advantage of being logically very simple. 



Appendix: Differential constraints 

We present two mathematical relations — equations ()A.1|) and (jA.4|) below — that have 
not been used in the body of the article, but that nevertheless are somewhat interest- 
ing by their own, since they represent differential constraints on momentum and on the 
Hamiltonian. 

Replacing (|3.6|) into ()3.5|) . one arrives at a single equation for p(u): 

dp k (u) 

Pi(u) = <pij{u) — u k . (A.l) 

d Uj 

On defining ipij(u) as the components of the inverse matrix of the (pij(u), that is ipij (fjk = 
5ik, equation (jA.lj) can be rewritten as 

^^ = ^WPr (A.2) 

Taking now the second derivative of Hamilton's equations (|3.9|) we get 

dm d 2 H 



dpj dpi dpj 



(A.3) 



Replacing (jA.3|) into equation (jA.2|) . we arrive at the following differential constraint on 
Hip): 

dH d 2 H 

7T = 7T-JT- V>;fc(%# p k . A.4 
dpi dpi dpj 
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